STABILITY OF RADIATIVE SHOCK PROFILES FOR 
HYPERBOLIC-ELLIPTIC COUPLED SYSTEMS 



TOAN NGUYEN, RAMON G. PLAZA, AND KEVIN ZUMBRUN 

Abstract. Extending previous work with Lattanzio and Mascia on the scalar (in fluid-dynamical 
variables) Hamer model for a radiative gas, we show nonlinear orbital asymptotic stability of small- 
amplitude shock profiles of general systems of coupled hyperbolic-eliptic equations of the type 
modeling a radiative gas, that is, systems of conservation laws coupled with an elliptic equation 
for the radiation flux, including in particular the standard Euler— Poisson model for a radiating 
gas. The method is based on the derivation of pointwise Green function bounds and description 
of the linearized solution operator, with the main difficulty being the construction of the resolvent 
kernel in the case of an eigenvalue system of equations of degenerate type. Nonlinear stability then 
follows in standard fashion by linear estimates derived from these pointwise bounds, combined with 
nonlinear-damping type energy estimates. 



1. Introduction 

In the theory of non-equilibrium radiative hydrodynamics, it is often assumed that an inviscid 
compressible fluid interacts with radiation through energy exchanges. One widely accepted model 
[37j considers the one dimensional Euler system of equations coupled with an elliptic equation for 
the radiative energy, or Euler-Poisson equation. With this system in mind, this paper considers 
general hyperbolic- elliptic coupled systems of the form, 

Ut + f(u) x + Lq x = 0, 

-q xx + q + g(u) x = 0, 

with (x, t) £ Rx[0, +oo) denoting space and time, respectively, and where the unknowns u £ U C R™, 
n > 1, play the role of state variables, whereas q € R represents a general heat flux. In addition, 
L 6 R™ xl is a constant vector, and / 6 C 2 {U\W a ) and g £ C 2 (U;M.) are nonlinear vector- and 
scalar- valued flux functions, respectively. 

The study of general systems like (jl.ip has been the subject of active research in recent years 
[lOl [TTJ [13l [17]. There exist, however, more complete results regarding the simplified model of a 
radiating gas, also known as the Hamer model [6], consisting of a scalar velocity equation (usually 
endowed with a Burgers' flux function which approximates the Euler system), coupled with a scalar 
elliptic equation for the heat flux. Following the authors' concurrent analysis with Lattanzio and 
Mascia of the reduced scalar model [16] , this work studies the asymptotic stability of general radiative 
shock profiles, which are traveling wave solutions to system (|1.1[) of the form 

u(x, t) = U(x - st), q(x, t) = Q(x - st), (1.2) 
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with asymptotic limits 

U(x) — » u±, Q{x) — > 0, as x — > ±00, 

being u± € U C M" constant states and set the shock speed. The main assumption is that the 
triple (m-|_,u_,s) constitutes a shock front |19j for the underlying "inviscid" system of conservation 
laws 

u t + f(u) x =0, (1.3) 
satisfying canonical jump conditions of Rankinc-Hugoniot type, 

/(u+) - /(u_) -«(«+- u_) = 0, (1.4) 

plus classical Lax entropy conditions. In the sequel we denote the jacobians of the nonlinear flux 
functions as 

A(u) := Df(u) e R" x ", B(n):=J} 9 (!i)el lx ", «6«. 

Right and left eigenvectors of A will be denoted as r € R nxl and I £ R lx ™, and we suppose that 
system (|1.3[) is hyperbolic, so that A has real eigenvalues a\ < • • • < a n . 

It is assumed that system represents some sort of regularization of the inviscid system (|1.3|) 
in the following sense. Formally, if we eliminate the q variable, then we end up with a system of 
form 

Ut + f(u) x = (LB(u)u x ) x + (u t + f(u) x ) xx , 

which requires a nondegeneracy hypothesis 

l p ■ (B ® L T r p ) > 0, (1.5) 

for some 1 < p < n, in order to provide a good dissipation term along the p-th characteristic field 
in its Chapman-Enskog expansion [34] . 

More precisely, we make the following structural assumptions: 

f,geC 2 (regularity), (SO) 

(SI) 



For all u 6 U there exists Aq symmetric, positive definite such that AqA 
is symmetric, and AqLB is symmetric, positive semi-definite of rank one 
(symmetric dissipativity =^ non-strict hyperbolicity) . Moreover, we assume 
that the principal eigenvalue a p of A is simple. 



No eigenvector of A lies in kei LB (genuine coupling). (S2) 



Remark 1.1. Assumption (|S1|) assures non-strict hyperbolicty of the system, with simple principal 
characteristic field. Notice that (|S1[) also implies that (Aq) 1 / 2 A(Aq)^ 1 / 2 is symmetric, with real 
and semi-simple spectrum, and that, likewise, (Aq) 1 / 2 B(Aq)^ 1 ^ 2 preserves symmetric positive semi- 
definiteness with rank one. Assumption (|S2[) defines a general class of hyperbolic-elliptic equations 
analogous to the class defined by Kawashima and Shizuta HH [35] and compatible with (|1.5[) . 
Moreover, there is an equivalent condition to (|S2[) given by the following 



Lemma 1.2 (Shizuta-Kawashima [TH [36]). Under (|S0[) - (|S1[) . assumption (|S2[) is equivalent to 
the existence of a skew- symmetric matrix valued function K :U — > R nx " such that 

Re (KA + A LB) > 0, (1.6) 

for all u £W. 

Proof. See, e.g., 0. □ 
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As usual, we can reduce the problem to the analysis of a stationary profile with s = 0, by 
introducing a convenient change of variable and relabeling the flux function / accordingly. Therefore, 
we end up with a stationary solution (U, Q) (x) of the system 

f(U) x + LQ X = 0, 
-Qxx + Q + g{U) x = 0. (L7) 

After such normalizations and under (|S0[) - (|S2[) . we make the following assumptions about the 
shock: 

/(u + ) = (Rankinc-Hugoniot jump conditions), (HO) 

P \ \ P (Lax entropy conditions), (HI) 

a p -i(u-) < < a p {u-), 

(Va p ) T r p ^ 0, for all u 6 U, (genuine nonlinearity), (H2) 



l P (u±)LB(u±)r p (u±) > 0, (positive diffusion). (H3) 

Remark 1.3. Systems of form arise in the study of radiative hydrodynamics, for which the 
paradigmatic system has the form 

Pt + {pu) x = 0, 
(pu) t + {pu 2 +p) x = 0, 

(p(e + 5« 2 )) t + (pu{e + \u 2 ) + pu + qj = 0, 

-q X x + aq + b(9 4 ) x = 0, 

which corresponds to the one dimensional Euler system coupled with an elliptic equation describing 
radiations in a stationary diffusion regime. In (|1.8|) . u is the velocity of the fluid, p is the mass 
density and 8 denotes the temperature. Likewise, p = p(p, 9) is the pressure and e = e(p, 9) is the 
internal energy. Both p and e are assumed to be smooth functions of p > 0, 9 > satisfying 

p p > 0, p e ^ 0, e e > 0. 

Finally, q = p\ x is the radiative heat flux, where x represents the radiative energy, and a, b > 
are positive constants related to absorption. System (|1.8|) can be (formally) derived from a more 
complete system involving a kinetic equation for the specific intensity of radiation. For this derivation 
and further physical considerations on (|1.8[) the reader is referred to [371 E] ■ 

The existence and regularity of traveling wave type solutions of under hypotheses (|S0[) - (|S2[) . 
(|H0p - (IH3|) is known, even in the more general case of non-convex velocity fluxes (assumption (|H2[) 
does not hold). For details of existence, as well as further properties of the profiles such as mono- 
tonicity and regularity under small-amplitude assumption (features which will be used throughout 
the analysis), the reader is referred to [T71 118j . 

1.1. Main results. In the spirit of [S1I22. 24, 25j, we first consider the linearized equations of (jl.lj) 
about the profile (U, Q): 

u t + (A(U)u) x + Lq x = 0, 

-q xx + q+(B{U)u) x = 0, [i - J > 
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with initial data u(0) — u . Hence, the Laplace transform applied to system (|1.9p gives 

Xu + (A(U) u) x + Lq x = S, 

~q xx + q+(B(U)u) x =0, (L10) 

where source S is the initial data uq. 

As it is customary in related nonlinear wave stability analyses (see, e.g., [U [33l ETJ [38] ) , we start 
by studying the underlying spectral problem, namely, the homogeneous version of system (|1.10[) : 

Xu + (A(U) u) x + Lq x = 0, 

-q xx +q+{B(U)u) x =Q. (L11) 

An evident necessary condition for orbital stability is the absence of L 2 solutions to for 
values of A in {Re A > 0}\{0}, being A — the eigenvalue associated to translation invariance. This 
spectral stability condition can be expressed in terms of the Evans function, an analytic function 
playing a role for differential operators analogous to that played by the characteristic polynomial for 
finite-dimensional operators (see [U [33J, O 2H [22] and the references therein) . The main property 
of the Evans function is that, on the resolvent set of a certain operator C, its zeroes coincide in 
both location and multiplicity with the eigenvalues of C. Thence, we express the spectral stability 
condition as follows: 

There exists no zero of the Evans function D on {Re A > 0} \ {0}; equiva- 
lently, there exist no nonzero eigenvalues of C with Re A > 0. 

Like in previous analyses 41. 38, 40 , we define the following stability condition (or Evans function 
condition) as follows: 



There exists precisely one zero (necessarily at A = 0; see Lemmas 12. 51 - [2^6]) ,„x 
of the Evans function on the nonstable half plane {Re A > 0}, 



which implies the spectral stability condition (|SS[) plus the condition that D vanishes at A = at 
order one. Notice that just like in the scalar case [16], due to the degenerate nature of system ([1.11)1 
(observe that A(U) vanishes at x — 0) the number of decaying modes at ±oo, spanning possible 
eigenfunctions, depends on the region of space around the singularity. Therefore the definition of 
D is given in terms of the Evans functions D± in regions x ^ 0, with same regularity and spectral 
properties (see its definition in (|2 . 23[) and Lemmas 12.51 - [2T61 below V 

Our main result is then as follows. 

Theorem 1.4. Assuming ([1.5j) . (|S0 |) - ([S2[) . (|H0[) - (|H3|) . and the spectral stability condition ([D]), then 
the Lax radiative shock profile (U, Q) with sufficiently small amplitude is asymptotically orbitally 
stable. More precisely, the solution (u, q) of (|1.1|> with initial data uq satisfies 

\u(x,t)-U(x-a(t))\ LP <C(l+t)"' (1_1/p) l"o|Lin^ _ 
\u{x,t) -U{x-a(t))\ H i < C(l + t)- 1/4 K|Linff* 

and 

\q(x,t) - Q(x - a(t))\ m * <C(l+*r l(1 - 1/p) kU lr w4 

\q(x,t) - Q(x - a(t))\ H * < C(l + ir 1/4 KUi nff 4 

for initial perturbation uq := uq — U that are sufficiently small in L 1 n H 4 , for all p > 2, for some 
a(t) satisfying a(0) = and 

\a(t)\ < C\u Q \ L i nH 4 

i/o (1-14) 
d(t|<Cl + t- 1/a tioUw- 
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Remark 1.5. The time-decay rate of q is not optimal. In fact, it can be improved as we observe 
that \q(t)\ L 2 < C\u x (t)\ L 2 and \u x (t)\ L 2 is expected to decay like i -1 / 2 ; however, we omit the 
detail of carrying this out. Likewise, assuming in addition a small L 1 first moment on the initial 
perturbation, we could obtain by the approach of [32] the sharpened bounds |d| < C(l + t) a ~ x , 
and la — a (+oo)| < C(l +ty- 1 / 2 , for a > arbitrary, including in particular the information that 
a converges to a specific limit (phase-asymptotic orbital stability); however, we omit this again in 
favor of simplicity. 

We shall prove the following result in the appendix, verifying Evans condition (D). 
Theorem 1.6. For e := \u+ — U—\ sufficiently small, radiative shock profiles are spectrally stable. 
Corollary 1.7. The condition (D) is satisfied for small amplitudes. 

Proof. In Lemmas 12.51 - [2TB1 below, we show that D{\) has a single zero at A = 0. Together with 
Theorem ll.6l this gives the result. □ 

1.1.1. Discussion. Prior to [16] , asymptotic stability of radiative shock profiles has been studied in 
the scalar case in [12] for the particular case of Burgers velocity flux and for linear g(u) = Mu, with 
constant M. Another scalar result is the partial analysis of Serre [35] for the exact Rosenau model. 
In the case of systems, we mention the stability result of [21] for the full Euler radiating system under 
special zero-mass perturbations, based on an adaptation of the classical energy method of Goodman- 
Matsumura-Nishihara 4 , 27] . Here, we recover for systems, under general (not necessarily zero-mass) 
perturbations, the sharp rates of decay established in [12] for the scalar case. 

We mention that works |12[ 116] in the scalar case concerned also large- amplitude shock profiles 
(under the Evans condition (D), automatically satisfied in the Burgers case [12]). At the expense 
of further effort book-keeping- specifically in the resolution of flow near the singular point and 
construction of the resolvent- we could obtain by our methods a large-amplitude result similar to 
that of [16] . However, we greatly simplify the exposition by the small-amplitude assumption allowing 
us to approximately diagonalize before carrying out these steps. As the existence theory is only for 
small-amplitude shocks, with upper bounds on the amplitudes for which existence holds, known to 
occur, and since the domain of our hypotheses in [16] does not cover the whole domain of existence 
in the scalar case (in contrast to [T2], which does address the entire domain of existence), we have 
chosen here for clarity to restrict to the small-amplitude setting. It would be interesting to carry 
out a large-amplitude analysis valid on the whole domain of existence in the system case. 

1.2. Abstract framework. Before beginning the analysis, we orient ourselves with a few simple 
observations framing the problem in a more standard way. Consider now the inhomogcncous version 

u t + (A(U) u) x + Lq x = g, 

-q xx + q+(B(U)u) x =h, ( ' ' 

of (|1.9p . with initial data u(x,0) = uo- Defining the compact operator K. := (—d x + 1) _1 of order 
— 1, and the bounded operator 

J := d x LKd x B{U) 
of order 0, we may rewrite this as a nonlocal equation 

u t + (A(U) u) x +Ju = d x LJCh + g, 

f n\ f \ ( L16 ) 
u{x, 0) = uq(x) 

in u alone, recovering q by 

q = -Kd x B{U)u + Kh. (1.17) 
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The generator C := —(A(U) u) x —Ju of (|1.16p is a zero-order perturbation of the generator —A(U)u x 
of a hyperbolic equation, so generates a C° semigroup e ct and an associated Green distribution 
G(x, t; y) := e S y (x). Moreover, e ct and G may be expressed through the inverse Laplace transform 
formulae 

e Ct = ±-.\ e xt (\ - £)~ 1 d\, 

Jr >-™° (1.18) 



G(x,t;y) = — e xt G x (x,y)d\ 
Jn—ioD 



'V 

for all r\ > 770, where G\(x, y) := (A — C)~ 1 5 y (x) is the resolvent kernel of C. 

Collecting information, we may write the solution of (|1.15[) using Duhamel's principle/ variation 
of constants as 

r+oa 

u(x,t) = / G(x,t;y)u (y)dy 



+ 00 



G(x, t - s; y){d x LK.h + g)(y, s) dyds, ( L19 ) 



q(x,t) = ((-JCd x B(U))u + JCh\(x,t), 

where G is determined through (|1.18[) . 

That is, the solution of the linearized problem reduces to finding the Green kernel for the it- 
equation alone, which in turn amounts to solving the resolvent equation for C with delta- function 
data, or, equivalently, solving the differential equation (|1.10[) with source 5* = S y (x). This we shall do 
in standard fashion by writing (|1.10[) as a first-order system and solving appropriate jump conditions 
at y obtained by the requirement that G\ be a distributional solution of the resolvent equations. 

This procedure is greatly complicated by the circumstance that the resulting (n + 2) x (n + 2) 
first-order system 

&{x, \)W X — A(x, X)W (1.20) 
is singular at the special point where A(U) vanishes, with O dropping to rank n + 1. However, 
in the end we find as usual that G\ is uniquely determined by these criteria, not only for the 
values ReA > ?7o > guaranteed by C°-semigroup theory/energy estimates, but, as in the usual 
nonsingular case [7] , on the set of consistent splitting for the first-order system (|1.20p , which includes 
all of {ReA > 0} \ {0}. This has the implication that the essential spectrum of C is confined to 
{ReA < 0} U {0}. 

Remark 1.8. The fact (obtained by energy-based resolvent estimates) that C — A is coercive for 
Re A > 770 shows by elliptic theory that the resolvent is well-defined and unique in class of distribu- 
tions for Re A large, and thus the resolvent kernel may be determined by the usual construction using 
appropriate jump conditions. That is, from standard considerations, we see that the construction 
must work, despite the apparent wrong dimensions of decaying manifolds (which happens for any 
ReA > 0). 

To deal with the singularity of the first-order system is the most delicate and novel part of 
the present analysis. It is our hope that the methods we use here may be of use also in other 
situations where the resolvent equation becomes singular, for example in the closely related situation 
of relaxation systems discussed in [22j [25] . 



2. Construction of the resolvent kernel 

2.1. Outline. In what follows we shall denote ' = d x for simplicity; we also write A{x) — A(U) 
and B{x) = B(U). Let us now construct the resolvent kernel for C, or equivalently, the solution 
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of (jl.lOp with delta- function source in the u component. The novelty in the present case is the 
extension of this standard method to a situation in which the spectral problem can only be written 
as a degenerate first order ODE. Unlike the real viscosity and relaxation cases [22j [23J [24l [25] (where 
the operator L, although degenerate, yields a non-degenerate first order ODE in an appropriate 
reduced space), here we deal with a system of form 

ew' = ao, \)w, 

where 

'A 

h, 



e 



is degenerate at x = 0. 

To construct the resolvent kernel we solve 

{Qd x -A(x,X))g x (x,y) = 6 y (x), (2.1) 

in the distributional sense, so that 

{ed x -A{x,\))g x (x,y)=0, (2.2) 

in the distributional sense for all x ^ y with appropriate jump conditions (to be determined) at 
x = y. The first entry of the three- vector Q\ is the resolvent kernel G\ of C that we seek. 

Namely Q\ , is the solution in the sense of distribution of system (jl.lOp (written in conservation 
form): 

(An)' = - (A + LB)u + Lp + S y {x) 

q' = Bu-p (2.3) 
p' = -<?• 

2.2. Asymptotic behavior. First, we study at the asymptotic behavior of solutions to the spectral 
system 

(A(x)u)' = -(A + LB{x))u + Lp, 

q' = B(x)u-p, (2.4) 

p = -q, 

away from the singularity at x = 0, and for values of A ^ 0, Re A > 0. We pay special attention to 
the small frequency regime, A ~ 0. First, we diagonalize A as 

(A- \ 

A := L p AR p = a p (2.5) 

\0 At) 

where A\ < — 9 < 0, A\ > 9 > 0, and a p £ K, satisfying a p (+oo) < < a p (— oo). Here, L p , R p are 
bounded matrices and L p R p = I. Defining v := L p u, we rewrite (|2.4I) as 

(A(x)vY = -(A + LB + L' p AR p )v + Lp, 

q' = Bv-p, (2.6) 

p' = -q, 

where L :— L p L and B :— BR p . Denote the limits of the coefficient as 

A± := lim A(x), B± := lim B(x)R p . (2.7) 

X — >±OC X — >±oc 

The asymptotic system thus can be written as 

W' = A ± (X)W, (2.8) 
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where W = (v,q,p) T , and 

A±(A) = I B± -1 I . (2.9) 




To determine the dimensions of the stable/unstable eigenspaces, let A G R + and A — > 0, +oo, 
respectively. The 2x2 lower right-corner matrix clearly gives one strictly positive and one strictly 
negative eigenvalues (this later will give one fast-decaying and one fast-growing modes). In the 
"slow" system (as |A| — » 0), eigenvalues are 

fJ ,f(X) = -X/af + 0(X 2 ), (2.10) 

where are eigenvalues of A± = A(±oo). Thus, we readily conclude that at x = +oo, there are 
p + 1 unstable eigenvalues and n — p + 1 stable eigenvalues. The stable S + (X) and unstable U + (X) 
manifolds (solutions which decay, respectively, grow at +oo) have, thus, dimensions 

dimZ7+(A) = p+ 1, 

+ (2-11) 
dim S (A) = n — p + 1, 

in Re A > 0. Likewise, there exist n — p+ 1 unstable eigenvalues and p stable eigenvalues so that the 
stable (solutions which grow at — oo) and unstable (solutions which decay at — oo) manifolds have 
dimensions 

dim U~ (A) = p, 

(2.12) 

dimS* (A) = n — p + 2. 

Remark 2.1. Notice that, unlike customary situations in the Evans function literature [Tl |4T| l3"l 
[221 [23l [33] , here the dimensions of the stable (resp. unstable) manifolds S + and S~ (resp. U + and 
U~) do not agree. Under these considerations, we look at the dispersion relation 

fr±(*C) = -if - A ± l (A + LB±)£, 2 A ± l = 0. 

For each £ g R, the A-roots of the last equation define algebraic curves 

xfio g <r(i + LB^yH-e + iA±t(i + e 2 )), £ G M, 

touching the origin at ^ = 0. Denote A as the open connected subset of C bounded on the left by 
the rightmost envelope of the curves A^(£), (6l. Note that the set {Re A > 0, A ^ 0} is properly 
contained in A. By connectedness the dimensions of U ± (X) and 5 ± (A) do not change in A G A. We 
define A as the set of (not so) consistent splitting pQ, in which the matrices A±(A) remain hyperbolic, 
with not necessarily agreeing dimensions of stable (resp. unstable) manifolds. 



Lemma 2.2. For each A G A, the spectral system (12. 8|) associated to the limiting, constant coeffi- 
cients asymptotic behavior of (|2.4p . has a basis of solutions 

e ^W*V±(A), .r. O.j ! n ■ ■>. 

Moreover, for |A| ~ 0, we can find analytic representations for fij and Vj , which consist of2n slow 
modes 

M ±(A) = -X/ af + 0(A 2 ), j = 2, n + 1, 

and four fast modes, 

A*f (A) = ±0f + 0(A), 
M± +2 (A) = T^f +2 + 0(A). 
where 9f and 6 n+2 are positive constants. 
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In view of the structure of the asymptotic systems, we are able to conclude that for each initial 
condition xq > 0, the solutions to (|2.4[) in x > xq are spanned by decaying/growing modes 

(2 . 13 , 

* + : = W',t P+2 ,-,«,t«). 



as x — > +00, whereas for each initial condition xq < 0, the solutions to (|2.4|) are spanned in x < xq 
by growing/decaying modes 

* : = {</ > n-p+3'->n+2}> 

as x — > —00. 

We rely on the conjugation lemma of [29] to link such modes to those of the limiting constant 
coefficient system (|2.8|) . 

Lemma 2.3. For \X\ sufficiently small, there exist growing and decaying solutions ip^x, A), <fif(x, X), 
in x 0, 0/ c/ass C 1 m a; and analytic in X, satisfying 

i,f{x,X)=e»t^ x V±{\){I + 0{e-^% 

, (2-15) 
</>±(x, A) = e"i W*Vf(A)(/ + 0(e-"l*l)), 

where < rj is the decay rate of the traveling wave, and fij and Vp are as in Lemma \2.'2\ above. 

Proof. This a direct application of the conjugation lemma of [29] (see also the related gap lemma in 

dim 121113]). □ 



2.3. Solutions near x ~ 0. Our goal now is to analyze system (|2.4p close to the singularity x = 0. 
To fix ideas, let us again stick to the case x > 0, the case a; < being equivalent. We introduce a 
"stretched" variable £ as follows: 

so that £(1) = 0, and £ — > +00 as x — » + . Under this change of variables we get 

. du I du 1 

u = — — 



dx a p (x) d£ a p (x) ' 

and denoting ' = d/d£. In the stretched variables, making some further changes of variables if 
necessary, the system (|2.6|) becomes a block-diagonalized system at leading order of the form 

z= ^ +o„(oe(0^ (2.16) 

where 0(£) is some bounded matrix and a is the (p,p) entry of the matrix A + LB + L' p AR p + A' , 
noting that 

a(0 >S >0, 

for some So and any £ sufficiently large or x sufficiently near zero. 

The blocks —al and are clearly spectrally separated and the error is of order 0(|a p (£)|) — > 
as £ — » +00. By the pointwise reduction lemma (see Lemma IB. II and Remark IB. 21 below), we can 
separate the flow into slow and fast coordinates. Indeed, after proper transformations we separate 
the flows on the reduced manifolds of form 

Z 1 = -aZ 1 +0(a p )Z u (2.17) 

Z 2 = 0(a p )Z 2 . (2.18) 
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Since —a < — 5 < for A ~ and £ > 1/e, with e > sufficiently small, and since a p (£) — > as 
£ — > +oo, the Zi mode decay to zero as £ — > +oo, in view of 

e -/«a(*)d* < e -(HeA+|«o)€ < 

These fast decaying modes correspond to fast decaying to zero solutions when x — > + in the 
original u- variable. The Z2 modes comprise slow dynamics of the flow as x — » + . 

Proposition 2.4. There exists < eo 1 sufficiently small, such that, in the small frequency 
regime A ~ 0, f/ie solutions to the spectral system (|2.4p m (— eo, 0) U (0, eo) are spanned by fast modes 




±e o ^x^0, (2.19) 
decaying to zero as x — > ± , and slowly varying modes 

±e ^x^ 0, (2.20) 




loj'i/l bounded limits as x — > . 

Moreover, the fast modes (|2. 19p decay as 



«L - - ( 2 -21) 




OdxrV^a;)) - 0, jVp, (2.22) 

asi-t 0^; /iere, ao is some positive constant and Uk p — {uk p i, ■■■,Uk p pi ■■■,Uk p n) T ■ 
2.4. Two Evans functions. We first define the following related Evans functions 

£>±(»,A):=dst(*+Wg *-)(», A), fory^O, (2.23) 

where $ ± are defined as in ([2~T3"|) , ([2~T4]) . and = (u± ,Q% ,pf V are defined as in f2~T^|) . Note 
that A: p here is always fixed and equals to n — p + 2. 
We first observe the following simple properties of D± . 

Lemma 2.5. For A sufficiently small, we have 

D ± (y,X) = (det^-Vi^AA + OQAl 2 ), (2.24) 

where 

A:=dct(r+ ■•• r+_j rf p+1 •■• r~ +1 -[u]j 

fqf qf\ (2-25) 

MJ^ft [ti] = u+ — it_ and eigenvectors of (A±)~ 1 (LB)±, spanning the stable/unstable subspaces at 
±00, respectively. 
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Proof. By our choice, at A = 0, we can take 

4>t(x,0)=<f>- +2 (x,0) = W x (x) (2.26) 
where W is the shock profile. By Leibnitz' rule and using (|2.26p . we compute 

d x D_(y, 0) = det (d x <f>+, 0+ _ x> W+, fa p+1 , </>- +2 ) + ■ ■ ■ 

• • • + dot (4>+, 0+ _ v W+,^ p+1 , dx^~ +a ) , 

where, by using (|2.26[> . only the first and third terms are possibly nonvanishing and thus grouped 
together, yielding 



x D_(y, 0) = det (fi, 0+ _ v W+,^ p+v d x ^ +2 - d x <f>+) ^ (2.27) 



Recall that Wu , cjyf satisfy 



BW X = A(z, X)W, (2.28) 

where W = (u,q,p) and 

' A 

h 



e 

Thus, 9a0^(x, A) satisfies 



&(d X 4>t) x = A(i, O)5 A 0+(:r, 0) + A A(ss, 0)</)+(x, 0), 

which directly gives 

(ad A uf ) x = -L(d x q+) x - Sx. (2.29) 
Likewise, d x <p~ +2 (x , A) = (d x u~ +2 ,d x q~ +2 ,d x p~ +2 ) satisfies 

(ad A u~ +2 )* = -i(9 A C+2)^ - ( 2 - 30 ) 

Integrating equations (|2.29|) and (|2.30|) from +oo and — oo, respectively, with use of boundary 
conditions d x <pi (+oo) = d x 4>~ +2 (—oo) = 0, we obtain 

Ad x ut = —Ld x qt — u + u + 

(2.31) 

Ad x u n+2 = -Ld x q n+2 -u + u— 

and thus 

Md X K +2 - d x u+) = -L{d x q- +2 - d x q+) - [u\. (2.32) 

In addition, we note that W£ , 4>f satisfy the equation (|2.28[) and thus (Ail)' = —Lq' with 
W fe +(+oo) = 4>t (+oo) = 0, 0~ +2 (-oo) = 0, qf (±oo) = 0, and 



uj (+oo) = (A + )- L rj, j = 2, ...,k p - 1 
«7(-oo) = (A-)- 1 ^, j = k p + l,...,n+l 
Thus, we integrate the equation (An)' = —Lq', yielding 

Au+ = -Lqf, for j = 1, ftp, 
Au+ = -Lq+ + r+, for j = 2, k p - 1, 



(2.33) 

AUj = —Lq- +r j , for j = k p + 1, n + 1 
y4uj = — Lg~, for j = n + 2. 
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Using estimates (|2.33p and (|2.32|) , we can now compute the A-derivative (|2.27p of D± at A = as 



tut 



d x D_(y,0) =det 



\Pi ■■■ Pj 

(deti4) _1 det 9i + 
W 



4 P 
pt 



r 



Pj 



q t 
Pt 



d\u~ +2 - d\u~l 



Pi 



d \Pn+2 ~ d Xpt 



(2.34) 



(det^) _1 det 



det ( r. 



K Pi ?t, 

which proves (|2.24|) . The proof for D + follows similarly. 
Lemma 2.6. Defining the Evans functions 

At (A) := At(±l,A), 

we then have 

D + (\)=mD-(\) + 0(\\\ 2 ) 

where m is some nonzero factor. 
Proof. Proposition 12.41 gives 



fep+i 



n+l 



W kA X ) = 



□ 

(2.35) 
(2.36) 

(2.37) 



as i — > 0, where olq is defined as in Proposition 12.41 which guarantees an existence of positive 
constants e\ , e 2 near zero such that 



Thus, this together with the fact that are solutions of the ODE (I2.28P yields 

">£(-!) = m ^ w k p (+ 1 ) 

for some nonzero constant Wfc . Putting these estimates into (|2.24p and using analyticity of D± in 
A near zero, we easily obtain the conclusion. □ 

3. Resolvent kernel bounds in low-frequency regions 

In this section, we shall derive pointwise bounds on the resolvent kernel G\(x, y) in low- frequency 
regimes, that is, |A| — > 0. For definiteness, throughout this section, we consider only the case y < 0. 
The case y > is completely analogous by symmetry. 

We solve (|2.3p with the jump conditions at x = y: 

(A(y)- 1 0\ 
[&(-,!/)]= 10 (3.1) 
V 1/ 

where, working on diagonalized coordinates (see (|2.6p ). we can assume that A is of diagonal form as 
in (1231). 



A, 



A 







A}, 
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with A 1 < —0 < 0, A%(y) > 8 > 0. Meanwhile, we can write G\(x, y) in terms of decaying solutions 
at ±oo as follows 

Qx ( x y)= h + (x,yC + (y,X) + W+(x,\)C+(y,\), x>y, 



where Cj are row vectors. We compute the coefficients by means of the transmission conditions 
(|3.ip at y. Therefore, solving by Cramer's rule the system 



fC+\ (Aiy)- 1 0\ 

(V W+ $-) \C+\ = 10, (3.3) 



we readily obtain, 



C+ j ( 2/ ,A)= J D_(y,A)- 1 (<i> + W+ f 1 oj (3.4) 

where M a * denotes the adjugate matrix of a matrix M. Note that 

Cf p (y,A)=a p (y)- 1 J D_( 2 / ) A)- 1 ($ + W+ A), (3.5) 

C±(y ) A) = ^D_( 2 /,A)- 1 ($+ W£ $-) fcJ (j/,A)(^)- 1 ) H , ^p, (3.6) 



where () lJ is the determinant of the minor, and (A(y) )klt I ^ P, are bounded in y. 
We then easily obtain the following. 

Lemma 3.1. For y near zero, we have 

ct(y,x)= i« o (M) + 0(i), 



C- +2 (y,\) = - T v ([u]) + O(l), 



1 
A 

where Vq([u]) is some constant vector depending only on [u] and 



(3.7) 



C+(y,\)=a p (y)- 1 \y\- a °0(l), 

C+(y,X)=0(l) Kj<k p , (3.8) 
Cr(y,\) =0(1) k p <j<n + 2, 

where k p = n — p + 2, c<q is defined as in Proposition^^ and 0(1) is a uniformly bounded function, 
probably depending on y and A. 

Proof. We shall first estimate C~ +2 (y, A). Observe that 



($ + W+ *-) (y,A) = (*+ W fc + (y,0) + O(A) 

where by the same way as done in Lemma 2.5 we obtain an estimate 

W+ *-) (y,0) = a p (detA)- 1 7 _(y)A^+ 2 , 
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where 7_(y) and A are denned as in (|2.25[) . and A p,n+2 denotes the minor determinant. Thus, 
recalling (|2.24[) and (|3.5p . we can estimate C~ +2p (y, A) as 

/ \ p,n+2 

C- +2 , p (y 1 \) = a p {yy 1 D-{vA)- l ^ + W+ *") (y, A) 

= -Ia-'a"'^ 2 + o(i), 

A 

where 0(1) is uniformly bounded since a p (y)D-(y, A) and normal modes <^~~ are all bounded uni- 
formly in y near zero. Similar computations can be done for C~ +2l {y, A). Thus, we obtain the 
bound for C~ +2 as claimed. The bound for follows similarly, noting that 4>~ +2 = 4>l at A = 0. 

For the estimate on , we first observe that by view of 1(2.240 . with noting that det(A) ~ a p (y) 
as \y\ — ► 0, and the estimate (|2.19[) on iv~£ p , 

\D_(y,\)\>6\\\\ y r, (3.9) 

for some 9 > 0. This together with the fact that </>~ +2 = 4>\ at A = yields the estimate for C£ as 
claimed. 

We next estimate (resp. C^~) for 1 < j < fe p (resp. k p < j < n + 2). We note that by view of 
estimate (|2.19|) on Wfe p , 

$+ $-) =0(A)0(|2/r o a p ( y )) 

and for k p, 

(V $-) feJ = 0(A)0(|yr) 

These estimates together with (|3 .9[) and (|3.6[l . (|3.5[) immediately yield estimates for as claimed. 

□ 

Proposition 3.2 (Resolvent kernel bounds as \y\ — > 0). For y near zero, there hold 

G x (x 1 y) = \- 1 W x v Q ([u]) + 0(l) e -^ a t + °^ x + 0(e- e ^) (3.10) 

for y < < x, and 

G x (x, y) = X- l W x v ([u]) + 0(1) (l + ^ ) (3.11) 

v a p ( y )\y\ a J 

for y < x < 0, and 

G x {x, y) = X-^MVA) + O(l) E e- (A/ ^ +0(A2))a; + ©(e-^l) (3.12) 

aT<0 

/or x < y < 0. 

Similar bounds can be obtained for the case y>0. 

Proof. For the case y < < x, using (|3.7|) and recalling that <fti{x) = W x + 0(A)e~ 9 ' x and 
Wjl (x) = 0, we have 

G x (x,y) = <P+(x)C+(y) = £ tf(x)C+(y) 

= (VF, + 0(A)e- e l^l) (^o(M) + 0(1)) + 0(1) £ e^, 

J=2 
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(3.13) 



yielding (|3.10|k here, we recall that 

l*f = -X/af + 0{\ 2 ) 

with > for j = 2, k p — 1 and a J < for j = k p + 1, n + 1 (a^ are necessarily eigenvalues 
of A± ). In the second case y < x < 0, from the formula (|3.2p . we have 

G x (x, y) = <P + (x, X)C+ (y, A) + W+ (x, \)C+ (y, A) 

where the first term contributes A _1 Wo(M)W x + 0(1) as in the first case, and the second term is 
estimated by $£8$ and (pT2"T) . 

Finally, we estimate the last case x < y < in a same way as done in the first case, noting that 
y is still near zero and W~ +2 (x) = W x + 0{X)e~ e ^ . □ 

Next, we estimate the kernel G\(x, y) for y away from zero. Note however that the representations 
(|3.2[) and above estimates fail to be useful in the y — > — oo limit, since we actually need precise decay 
rates in order to get an estimate of form 

\Gx{x,y)\<Ce-^-v\, 

which are unavailable from (j>f in the y — > — oo regime. Thus, we need to express the (+)-bases in 
terms of the growing modes ipj at — oo, and the decaying mode <p~ where ipj, cf>J are defined as in 
Lemma 12.31 Expressing such solutions in the basis for y < 0, away from zero, there exist analytic 
coefficients djk(X), ejk(X) such that 

0+ 0, A) = ^2 d Jk 0; A ) + e i k A ) 
W+{x,X) = ^4^(A)C(^ A) l^e M (%-(^,A). 
Furthermore, for our convenience, we define the following adjoint normal modes 

$-) _1 9" 1 . (3.14) 

We then obtain the following estimates. 

Lemma 3.3. For |A| sufficiently small and \x\ sufficiently large, 

$r{x, A) = 0(e-^W*)Vr(X)(I + O^^)), 

4>J(x, A) = 0(e-^^ x )Vr(X)(I + 0(e~ e ^)) 
where fij are defined as in Lemma \2.3[ 

Proof. The proof is clear from the estimates of ipj, <fij in (|2.15|) . □ 
Lemma 3.4. We have 

Cf(y,X)=J24k(Wk(y^T (3-16) 
Cj{y, X)=J2 cJ k {X)^(y, A)* + 4>J(y, A)*, (3.17) 
for meromorphic coefficients c^ k in A. 

Proof. The proof follows by using (|3.13[) . definition (|3.14[) . and property of computing determinants. 

□ 

We then have the following representation for G\(x,y), for y large. 



(3.15) 
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Proposition 3.5. Under the assumptions of Theorem \l-4\ for |A| sufficiently small and \y\ suffi- 
ciently large, we have 

G A (x,y) = ^ C +(A)0+( a; ,A)^7(y,A)*, (3.18) 
j.k 

for y < < x, and 

G x {x,y) = Y, d% (A)0T (a-, X )fa ( y , A)* - £ C W (y, A)*, (3.19) 

j,k k 

for y < x < 0, and 

j.k k 

for x < y < 0, where c^ k (X),d^ : k (X) are scalar meromorphic functions satisfying 

-l 



c + = (-I kp 0)(«+ W£ 

(o -i n -k P )(^ + w+ $-)~v. 



Proof. Using representation (|3.2p of G x (x,y) together with (|3.16p and (|3. 1T|) . we easily obtain the 
expansions (|3.18[) and (|3.20[) . respectively. For (|3.19[) , again, using (|3.16[) , (|3. 13[) . and (|3.2[) . we can 
write 

g x (x, y) = j2 d%{\)^(x, \)fa(y, \y+J2 e t^7( x > x r 

j,k j.k 



(3.21) 



Meanwhile, by (O and (|3l| . 

G A (x,y)=($+ W£ o) (*)(*+ W£ $-) _1 ( y )e- 1 (y) (3.22) 
In view of the definition (|3~T4|) of (f3~2Tj) and (f3~22|) yield 

(*- $-)e(i> + w+ o)(< 

' i - (o *-) (#+ w+ 



<J>-i <J> ) " $ 



o j U 



)( 

In-k p J V 



<J> + W+ <T> ) " 



which proves the proposition. □ 

Proposition 3.6 (Resolvent kernel bounds as \y \ — > +oo). Make the assumptions of Theorem \1.4\ 

Then, for \y\ large, defining 

E x (x,y) = \-i J2 V- e^+°^yW x (x), 

a7> 

there hold 

G x (x,y) = E x (x,y) 

+ 0(l)( J2 e^ a ^ + °^y + 0(e- e ^))( ]T e HKT+c(A 2 ))x + 0(e- e W)) (3-23) 
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for y < < x, and 

G x (x,y) =E x (x,y) + <D(l) £ eC-V-r+WX'-v) 

aj>0 

+ 0(1) V e (Va7+e>(A 2 ))y e (-A/a-+0(A 2 )):z + (g-fl(|a:-v|)) 
aT>0, <0 

for y < x < 0, and 

G x (x,y) =E x (x,y) + 0(l) £ e (-V«7+o(A 2 ))(^) 



(3.24) 



a- <0 



+ 0(1) V e (A/ °7 +C(A 2 ))y e (-A/a-+O(A 2 )) :c + c)( | e ~0(| a: - y |)^ 



(3.25) 



a, >0, a. <0 



/or x < y < 0. 

Similar bounds can be obtained for the case y>0. 

Proof. The proof follows directly from the representations of G x (x, y) derived in Proposition l3.5l and 
the corresponding estimates on normal modes, noting that 



\ c jk\'\ a jk\ I (l) otherwise. 



Indeed, we recall, for instance, that 

kj 



c+ k = DZ l {-h p 0)(V W+ $-) #" 



where () ■? denotes the determinant of the (k,j) minors. For the case j ^ 1, by using the fact that 
we choose <f>^ = <fi~ +2 = W x at A = 0, determinant of the (fc, j) minor therefore has the order one in 
A, which cancels out the A -1 term coming from our spectral stability condition: IDZ 1 ] < 0(A -1 ). 

□ 

4. POINTWISE BOUNDS AND LOW-FREQUENCY ESTIMATES 

In this section, using the previous pointwise bounds (Propositions 13.21 and 13. 6|) for the resolvent 
kernel in low-frequency regions, we derive pointwise bounds for the "low-frequency" Green function: 

G\x,t;y):=±- [ e xt G x (x,y)d\ (4.1) 

lm Jrn{|A|<r} 

where T is any contour near zero, but away from the essential spectrum. 



Proposition 4.1. Under the assumptions of Theorem \ l-4\ defining the effective diffusion (3± := 
(L p LBRp)± (see (|2.5|) ), the low-frequency Green distribution G 1 (x, t; y) associated with the linearized 
evolution equations may be decomposed as 

G I (x,t;y)=E + G I + R, (4.2) 

where, for y < 0: 

E(x,t;y):= ]T U s (x)V£ e k (y,t), (4.3) 

a->0 

■*•"=-(-* (^) (44) 
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\d2dl& {x,t;y)\ < Ct-M+hMfet-We 

k—l 

, V V f -i/2 p -(x-a-(t-\ V /a-\)) 2 /Mt\ 



? -{*-V-*l t) 2 /Mt 

' k=l 



(4.5) 



a. <0,a. >0 



(4.6) 



R(x, t; y) = 0(e-^*-yW) + O^Mx, y) 1 + —-(x/y) a 

1 a P {y) 

for some t), C, M > 0, where < \/3\, \j\ < 1, a = LB \^,^f 0) and 

X(xy) = { 1 > - 1 <y< x <° 
x ^' y > \0, otherwise. 

Symmetric bounds hold for y > 0. 

Proof. Having the resolvent kernel estimates in Propositions 13.21 and 13 - 61 we can now follow the 
previous analyses of [JTJ [521 H3] • Indeed, the claimed bound for E precisely comes from the A -1 
term. Likewise, estimates of G 1 are due to bounds in Proposition ^. 61 for y away from zero and those 
in Proposition 13.21 for y near zero but x away from zero. The singularity occurs only in the case 
— 1 < y < x < 0, as reported in Proposition 13.21 In this case, using the estimate (|3.1ip and moving 
the contour T in (|4.1|) into the stable half-plane {Re A < 0}, we have 



( i +J£>= o( ^ i+ J£^ 



a P {y)\y\ a ) a p {y)\y\«)' 

which precisely contributes to the second term in R(x,t\y). The first term in R(x,t;y) is as usual 
the fast decaying term. □ 

With the above pointwise estimates on the (low-frequency) Green function, we have the following 
from [HH3I. 

Lemma 4.2 ([22j[23]). Under the assumptions of Theorem \1.4[ G 1 satisfies 



d^G I (;t;y)f(y)dy <C(l+t)^ (1/9 ~ 1/p) ^ l/2 |/k,, (4.7) 



for all t > 0, some C > 0, for any 1 < q < p. 

We recall the following fact from [55] . 
Lemma 4.3 ( 39 ). The kernel e satisfies 

\e y {;t)\ LV ,\e t (;t)\ L „<Cr^ l -^\ 
\e yt (-,t)\L,<Cr^~^-^. 
for all t > 0, some C > 0, for any p > 1. 

Finally, we have the following estimate on R term. 
Lemma 4.4. Under the assumptions of Theorem \1.4\ R(x,t;y) satisfies 



(4.E 



R(;t;y)f(y)dy < C e -^(\f\ LP + |/| £ „.), (4.9) 

-or, LP 



for all t > 0, some C, r\ > 0, for any 1 < p < oo. 
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Proof. The estimate clearly holds for the fast decaying term e _ ''(l 2; -J'l+*) j n Whereas, to estimate 
the second term, first notice that it is only nonzero precisely when — l<y<:r<0or0<a;<y<l. 
Thus, for instance, when — 1 < x < 0, we estimate 



1 + 



f 



a P (y) 



<C\f\ L o 
< C\f\ L o 



(x/yT f(y)dy 



1 



1 



f 



a P (y) 



(x/y) a f(y)dy 



-l My) 



■{x/y) a dy 



where the last integral is bounded by that fact that a p (x) 
obtain 



,-vt 



x(x,y) 



i 



l 



a p{y) 



(x/y) a f(y)dy 



LP(-1,0) 



x as \x\ — > 0. From this, we easily 

<Ce-*\f\z-, 



which proves the lemma. 



□ 



Remark 4.5. We note here that the singular term a' p 1 (y)(x/y) a appearing in (|3.11j) and (14. 6j) 
contributes in the time-exponential decaying term. This thus agrees with the resolvent kernel for 
the scalar convected-damped equation ut + a p u x = —LBu, for which we can find explicitly the Green 
function as a convected time-exponential decaying delta function similar as in the relaxation or real 
viscosity case. 



5. Nonlinear damping estimate and high-frequency estimate 

In this section, we establish an auxiliary damping energy estimate. We first recall the nonlinear 
perturbation equations with (u, q) perturbation variables 

u t + (A(u)u) x + Lq x = a(U x + u x ), 

(5.1) 

-Qxx + q+ (B(u)u) x = 0, 

where we now denote 

A(u) :=Df(U + u), B(u) := Dg{U + u). (5.2) 

We prove the following: 



Proposition 5.1. Under the assumptions of Theorem \l-4\ s0 long as ||u||^2,oo and \a\ remain 
smaller than a small constant £ and the amplitude \U X \ is sufficiently small, there holds 

\\u\\ 2 Hk (t) < Ce- et \\u\\ 2 Hk (0) + C [ e- e( *- s) (IMIi 2 + \a\ 2 )(s)ds, 9 > 0, (5.3) 



for k= 1,...,4. 

Proof. Let us work for the case a = 0. The general case will be seen as a straightforward extension. 
We first observe that 

\A 0x \, \A ot \, \A X \, \A t \, \B X \, \B t \ = OQU x \ + Q (5.4) 



where A, B arc defined as in (15. 2j) and Aq the symmctrizer matrix as in ([SI 
We note that from the second equation of (|5.1|) we easily obtain 



\\q\\H* <C\\u\\ Hk -r, (5.5) 
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for k > 1. Meanwhile, from the first equation, we estimate 
\jfc( A ° u ' u ) = ( A o u t,u) + ^(A ot u,u) 

= -(A A x u + A Au x + Lq x ,u) + -(A ot U,u) 
= -(A A x u- -(A Q A) x u + Lq x ,u) + -(A Qt u,u), 



which, by (I5.4[) and (|5.5p . yields 



lj t (A u,u)<C\\u\\ 2 L2 . (5.6) 



Now, to obtain the estimates (|5.3p in the case of k = 1, we compute 

^(^0'"x,«a : ) = ((^0"i)x,W2;) + ~(A 0t U x ,U x ) - (A 0x U t ,U x ) 

= -{{A Q Au x + A Lq x ) x ,u x ) + ~(A ot u x ,u x ) - (A 0x u t ,u x ) 

= -(A Lq xx ,u x ) - (A Au xx ,u x },+(O(\U x \ + ()u x ,u x ) + || 
= -(A Lq xx ,u x ) + (O(\U x \+0u x ,u x )+O(l)\\u\\ 2 L2 

2 



l±(A d k u,d k u) = (A d k x u u d$u) + l -(A ot d k x u,d k x u) 



(5.7) 



(5.8) 



= -(A LBu xl u x ) + (0(\U X \ + CK, u x ) + 0(l)||u||| a , 
noting that since AqA is symmetric, we have 

-(A Au xx ,u x ) = -{(A A) x u x ,u x ) = (0(\U X \ + ()u x ,u x ). 
Likewise, in spirit of Kawashima-type estimates, we compute 

\~u( Ku ' Ux ) = \{ K t u -> u x) + ^(Ku t ,u x ) + ^(Ku,u xt ) 

= ^(K t u,u x ) + -(Ku t ,U x ) - ^(Ku x ,u t ) - ^(K x u,u t ) 
= (Ku t ,u x ) + ^(K t u,u x ) - ^(K x u 7 u t ) 

= - {KAu x + KA x u + KLq Xl u x ) + ^-(K t u,u x ) - ^(K x u,u t ) 

= -(KAu x ,u x ) + (<D(\U X \ + ()u x ,u x ) + 0(l)\\u\\ 2 L2 . 
Adding (|5.7|) and (|5.8|) together, we obtain 

~((Ku,u x ) + (A u x ,u x )) (59) 

= -((KA + A a LB)u x ,u x ) + (0(\U X \ + Qu x ,u x ) + C(l)||n||| 2 

which, by the Kawashima-type condition (|1.6|) : KA + AqLB > 9 and the fact that 0(\U X \ + £) is 
sufficiently small, yields 

1 d f,„ x . , , , \ 1. 

Similarly, for k > 1, paying attention to the leading terms, we can compute 



{{K u, u x ) + {A u Xl u x )^j < —-6{u x ,u x ) + (D(l)\\u\\ 2 L 2 (5.10) 



(d k x (A u t ),d k x u) + (0(\U X \ + C)d k x u,d k x u) + 0(l)\\u\\ 2 Hk ^, 
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where by using the first equation and then the second one, we obtain 

(d k (A u t ),d k u) = -(d k (A Au x + A A x u + A Lq x ),d k u) 

= -<4>£fi£ _1 fex, - (A Ad k+1 u, d k x u) +■■■ 

= -{A Q Ld k (Bu), d k x u) + ~((A A) x d*u, d k u) + ■■■ . 

Thus, we have obtained 

~(A d k u,d k u) , x 

2dV x ' x 1 (5.11) 

= -(A LBd k u,d k u) + (0(\U x \+()d k u,d k u) + 0(l)||u||^»_x. 
Meanwhile, we have the following fc t/l -order Kawashima-type energy estimate 

^{Kdt\ d k u) = (Kd k x - X u t , d k u) + ±(K t d k x -\ d k u) - l(K x dt l u, dt l u t ) 

= -{KAd k u,d k u) + (0(\U X \ +Od k u,d k x u) + 0(l)||«||| rt _ 1 . 

Hence, as before, adding (|5 . and (|5.12|) together and using the Kawashima-type condition 
((TTS|) : KA + A LB >9 and the fact that 0(\U x \ + () is sufficiently small, we obtain 

— ((Kdt^d^u) + (A d k u,d k u)) < - l -e(d k u,d k x u)+0{\)\\u\\ 2 Hk ^. (5.13) 
Now, for S > 0, let us define 

s 

£{t) -^{{Kdt^d^u) + (A d k u,d k xU y 

By applying the standard Cauchy's inequality on {Kd x x u, d k u) and using the positive definiteness 
of Aq, we observe that £ (t) ~ ||u||^ fe . We then use the above estimates (|5.6I) . (|5.10[) . (|5.13[) . and take 
5 sufficiently small to derive 

±£(t)<-e 3 £(t) + C\\u\\Ut) (5.14) 
for some #3 > 0, from which (|5.3[) follows by the standard Gronwall's inequality. □ 



With the damping nonlinear energy estimates in hands, we immediately obtain the following 
estimates for high-frequency part of the solution operator e ct : 

S 2 (t) = — x^mX^e^iX-C^dX, (5.15) 

for small positive numbers 0i,02] see (I1.18[) . Here, X{|imA|>e 2 } equals to 1 for |ImA| > Qi and zero 
otherwise. 

Proposition 5.2 (High-frequency estimate). Under the assumptions of Theorem \1.4\ 

||S 2 (f)/IU*<Ce-^||/|| ff2 , 
\\d^S 2 (t)f\\ L , <Ce-^\\f\\ Ha+2 , 

for some 0\ > 0. 
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Proof of the proposition follows exactly in a same way as done in our companion paper [16] for 
the scalar case. We recall it here for sake of completeness. The first step is to estimate the solution 
of the resolvent system 

Xu + (A u) x + Lq x = ip, 
-Qxx + q + (Bu) x = tp, 
where A(x) = Df(U) and B(x) = Dg{U) as before. 

Proposition 5.3 (High-frequency bounds). Under the assumptions of Theorem ] for some R,C 
sufficiently large and 7 > sufficiently small, we obtain 

I (A - C)-\y - Ld x (K,i/}))\ H i < c(\<p\ 2 H1 + |^ 2 ), 
| (A - £)~ x {<p - Ld x (Kip))\ L 2 < ^(Miri + M£ a ), 
for all |A| > R and Re A > -7, where K := {-d 2 x + l)" 1 . 



Proof. A Laplace transformed version of the nonlinear energy estimates (|5.3[) in Section[5]with k = 1 
(see [40], pp. 272-273, proof of Proposition 4.7 for further details) yields 

(ReA + y)|u|hi < c(l u li 2 + I^Ih 1 + IV"Il 2 )- (5-17) 

On the other hand, taking the imaginary part of the L 2 inner product of U against Xu = Cu + 
d x LICh + f and applying the Young's inequality, we also obtain the standard estimate 

|ImA||«|| a < \(Cu,u)\ + \(LKiji,u x )\ + \(<p,u)\ 

<c{\u\ 2 H1 +\^\ 2 L2 + \<p\%), 

noting the fact that £ is a bounded operator from H 1 to L 2 and K. is bounded from L 2 to H 1 . 
Therefore, taking 7 = 71 /4, we obtain from (|5.17|) and (|5. 18|) 



\MMh + W\m<c{\u\ 2 L2 + M 2 L , + \<p\%i), 



for any Re A > —7. Now take R sufficiently large such that \u\ 2 L2 on the right hand side of the above 
can be absorbed into the left hand side for |A| > R, thus yielding 

|A|M| 2 + M^ <c(m 2 l2 + M^), 

for some large C > 0, which gives the result as claimed. □ 
Next, we have the following 



Proposition 5.4 (Mid-frequency bounds). Under the assumptions of Theorem \1.4 

\(X - £y x <p\ L 2 < C\<p\ H i for R- 1 < \X\ < R and RcX > -7, 
for any R and C = C(R) sufficiently large and 7 = j(R) > sufficiently small. 



Proof. Immediate, by compactness of the set of frequency under consideration together with the 
fact that the resolvent (A — C)^ 1 is analytic with respect to H 1 in A; see, for instance, [39]. □ 



With Propositions 15.31 and 15.41 in hand, we are now ready to give: 
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Proof of Proposition \5. C A The proof starts with the following resolvent identity, using analyticity on 
the resolvent set p(C) of the resolvent (A — for all ip G £>(£), 

(A - = A _1 (A - C^Ccp + \- V 

Using this identity and (I5.15|) . we estimate 



Si{t)ip 



-71 -\-ioo 



71 —zoo 

71+ioo 



1 

1 

+ 2ri 
=: Si + S2, 

where, by Propositions 15.21 and 1 5 . 4\ we have 

f —71+200 

\Sx\ L 2 <C ' 



X { | ImA |> 72} e At A V^A 



| A |-l e RcA t | (A _ £) -l £ ^| i2 | dA | 

— 7l — ZOO 

/ — 7l 
\\\- 3/2 \Cv\ H i\d\\ 
-71—100 



and 



—71 +ZOO 



1 




+ 2^ 





-71+ir 
-71— ir 



L 2 



A'V^A 

-71 — 2.00 

<Ce-^\<p\ L *, 

by direct computations, noting that the integral in A in the first term is identically zero. This 
completes the proof of the bound for the term involving ip as stated in the proposition. The estimate 
involving t/j follows by observing that Ld x JC is bounded from H s to H". Derivative bounds can be 
obtained similarly. □ 



Remark 5.5. We note that in our treating the high-frequency terms by energy estimates (as also 
done in [151 1301 [It)] ), we are ignoring the pointwise contribution there, which would also be convected 
time-decaying delta functions. To see these features, a simple exercise is to do the Fourier transform 
of the equations about a constant state. 

6. Nonlinear analysis 

In this section, we shall prove the main nonlinear stability theorem. The proof follows exactly 
word by word as in the scalar case [16]. We present its sketch here for sake of completeness. Define 
the nonlinear perturbation 

U ](x,t):=(^\(x + a(t),t)-(^\(x), (6.1) 

where the shock location a(t) is to be determined later. 

Plugging (|6.ip into (|1.1|) . we obtain the perturbation equation 

u t + (Au) x + Lq x = Ni(u) x + a(t)(u x + U x ), 

(6.2) 

-q_xx + q + (Bu) x = N 2 {u) x , 

where Nj(u) = 0(|u| 2 ) so long as u stays uniformly bounded. 
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We recall the Green function decomposition 

G(x,t;y) = G I (x,t;y) + G II (x,t;y) 
where G I (x,t;y) is the low-frequency part. We further define as in Proposition 14. 11 

G T (x, t; y) = G 1 {x, t; y) - E(x, t; y) - R(x, t; y) 

and 

G u (x,t;y) = G u (x,t;y) + R(x,t;y). 
Then, we immediately obtain the following from Lemmas 14.21 14.41 and Proposition [52 



(6.3) 



Lemma 6.1. We obtain 
f +00 



t; y)f(y)dy < C(l + ^(^-i/p)-^/^^ 

LP 



for all 1 < q < p, [3 = 0, 1, and 



for all 2 < p < oo. 



G H (x,t;y)f(y)dy < Ce^f]^, 

Lp 



(6.4) 



(6.5) 



Proof. (|6.4| is precisely the estimate (|4. T[) in Lemma I4~2| recalled here for our convenience. (|6.5jl is a 
straightforward combination of Lemma l4.4l and Proposition l5.21 followed by a use of the interpolation 
inequality between L 2 and L°° and an application of the standard Sobolev imbedding. □ 

We next show that by Duhamel's principle we have: 
Lemma 6.2. We obtain the reduced integral representation: 



u(x, t) 



(G 1 + G H )(x, i; y)u (y)dy 

t p-\-oo 



J-oo 
t /* + oo 



!0 J-oo 

q(x,t) ={ICd x )(N 2 (u) - Bu){x,t) 



G*(x,t- s\y)(d y LK,N 2 {u) + Nx(u) + a(t)u) (y,s)dyds 
G n \x,t - s;y)(d y LK.N 2 {u) + Ni(u) + a(t)uj {y,s)dyds, 



(6.6) 



and 



a(t) 



d(i) 



-\-oo 

e t (y,t)u (y)dy 

— OO 

/ e y (y,t- s)(d v LJCN 2 (u) + Nx(u) + a(t)u)(y,s) dy ds. 

J-oo 
-\-oo 

e f (y,t)u (y)dy 

-OO 

t /— f-QO 



(6.7) 



e yt (y, t - s) {d y LKN2{u) + N^u) + a(t)uj [y, s) dy 



(6. 



ds. 
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Proof. By Duhamel's principle and the fact that 

G(x,t;y)U x {y)dy = e ct U x (x) = U x (x), 

J —OQ 

we obtain 

/+oo 
G(x,t;y)u (y)dy 
-oo 

ft />+oc 



J -oo 

a(t)U x . 



G(x, t - s; y) {d y LK,N 2 {u) + Ni (u) + a(t)uj (y, s) dy ds 



(6.9) 



Thus, by defining the instantaneous shock location: 

r+oo 

<*(*) = - / et(y,t)u (y)dy 



+oo 

e y (y, t — s) I dyLKNiyu) + N\{u) + a(t)u j (y, s) dy ds 



and using the Green function decomposition (|6.3[) . we easily obtain the integral representation as 
claimed in the lemma. □ 

With these preparations, we are now ready to prove the main theorem, following the standard 
stability analysis of [24 j [38 j [39] : 



Proof of Theorem \l-4\ Define 

(6.10) 



C(t):= sup \ u (s)\ LP (l + s )Hi-i/p) + l a ( s )\ + \a( s )\(l + s) 1 / i 

0<s<t,2<p<oo L 



We shall prove here that for alH > for which a solution exists with £(t) uniformly bounded by 
some fixed, sufficiently small constant, there holds 

at)<c(\u \ L i nH s+at) 2 ). (6.ii) 

This bound together with continuity of implies that 

C(*) < 2C\u \ L i nH s (6.12) 

for t > 0, provided that |tto|i,inir s < 1/4C 2 . This would complete the proof of the bounds as claimed 
in the theorem, and thus give the main theorem. 

By standard short-time theory/local well-posedness in H s , and the standard principle of continu- 
ation, there exists a solution u G H s on the open time-interval for which \u\h s remains bounded, and 
on this interval £(i) is well-defined and continuous. Now, let [0, T) be the maximal interval on which 
\u\h s remains strictly bounded by some fixed, sufficiently small constant S > 0. By Proposition [O] 
and the Sobolev embeding inequality |u|vi^2,oo < C\u\h°, s > 3, we have 



\u{t)\ 2 HS < Ce-"\uo\ a H . +C j\- e ^{\u(r)\l 2 + \af)dr 

<c(Ki 2 H3 + c(t) 2 )(i+tr 1/2 . 



(6.13) 



and so the solution continues so long as £ remains small, with bound (|6.12[) , yielding existence and 
the claimed bounds. 
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Thus, it remains to prove the claim (|6.1ip . First by representation (|6.6[) for u, for any 2 < p < oo, 
we obtain 



+00 



\u\ LP (t)< / (G 1 ^G II )(x 1 t;y)u Q (y)dy 



LP 



+ 00 



+ 00 



t —00 

--h+h + h, 



& y {x, t-s;y) (d y LK.N 2 (u) + Ni(u) + a(s)uj (y, s) dy 
G u (x, t-s;y) (d y LJCN 2 (u) + N^u) + d(t)u) (y, s) dy 



ds 



ds 



L" 



where estimates (|6.4[) and (|6.5[) yield 



LP 



/+00 ____ ____ 
(G 1 + G H )(x, t; y)u (y)dy 
-OO 

<Cr(l + t)-» (1 - 1/p) |«oUi + Ce-^|«o|H8 



and, with noting that d y LIC is bounded from L 2 to 1/ 



+00 



Gy(x, t-s-y) (d y L1CN 2 {u) + JVi(u) + d(t)«J (y, s) dy 

< C / (t - sJ-^-Vri-Vaflui^ + |d|)| u | La ( a )d8 

Jo 

< CCW 2 At - s )-i(l/2-l/p)-l/2 (1 + s )-3/4 ds 

Jo 

< CC(i) 2 (l + i)"5( 1 - 1 /p) 7 



ds 



and, together with (I6.13p . s > 4, 



/>+oo 



<C e 



< C 



G H {x, t-s;y) (d y LJCN 2 (u) + N^u) + a(s)uj (y, s) dy 
''( t_s ) \d y LKN 2 {u) + Ni(u) + a(t)u\ H 4 (s)ds 
^- s \\u\ H s + \a\)\u\ H s(s)d S 



ds 



LP 



<C(\u \ 2 Hs +((t) 2 ) / e -"(*- s )(l + s )- 1 d S 
Jo 

<c(\ Uo \%s+at) 2 )(i+t)- 1 . 



Thus, we have proved 



|«(t)Ml + < C(|«oUx n ^ + CW 2 )- 



(6.14) 



(6.15) 

Similarly, using representations (|6.7p and (|6.8[) and the estimates in Lemma 14.31 on the kernel 
e(y,t), we can estimate (see, e.g., [2~4l 155] ). 



|a(f)|(l + t) lj2 + \a(t)\ < C(\u \ L i + C(t) 2 ). 



(6.16) 
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This completes the proof of the claim (|6.11|) , and thus the result for u as claimed. To prove the 
result for q, we observe that K,d x is bounded from LP — > W 1,p for all 1 < p < oo, and thus from the 
representation (|6. 6|) for q, we estimate 

\q\w^(t) < C(\N 2 (u)\ LP + \u\ LP )(t) 

Ui 1/ n ( 6 - 17 ) 

< C\u\ LP (t) < C\u \ L , nH ,{\ + tJ-i^-Vp) 

and 

|9| HS+1 (f)<CH fls (t)<C|uo| £lnff »(l+i)- 1/4 , (6.18) 
which complete the proof of the main theorem. □ 

Appendix A. Spectral stability in the small-amplitude regime 

In this section we verify the spectral stability condition for small-amplitude profiles. Denoting 
A = A(U(x)), B = B(U(x)) we have the associated linearized spectral problem 

Xu + (Au) x + Lq x = 0, 

-q X x + q + (Bu) x = 0. 

Using the zero-mass conditions 

idx = 0, qdx = 0, 



we recast system (IA.1I) in terms of the integrated coordinates, which we denote, again, as u and q. 
The resulting system reads 

Xu + Au x + Lq x = 0, (A. 2) 

~q xx + q + Bu x = 0. (A. 3) 

In what follows we assume that the shocks are weak, that is, u± e A/"(w*), being Af a neighborhood 
of a certain state it*, for which 

< max \u — ttJ < e < 1, 

with e > sufficiently small; clearly, 

\u* — u±\, \u- — u + \ = 0(e) 
and the shock profile for U is approximately scalar, satisfying, 

U x = 0(e 2 )e-^(r p (u*) + 0(e)), 
U xx = 0(e 3 )e- 9 ^, <: ' 41 

for some 9, r\ > 0. For the principal characteristic field a p := a p (U(x)) we have 

( a p)x = 0(U X ) < 0, (monotonicity), (A. 5) 

We shall make use of the following 
Lemma A.l. Under (|S0|) - (|S2[) . there exists a scalar function (3 = f3(u) > 0, such that 

(A Q L) T = 0B, (A.6) 

for all u GW. 

Proof. Follows by elementary linear algebra facts, since AqLB is positive semi-definite with rank one 
and can be written as z <g> w, for some vectors z and w. It follows the existence of a scalar (3, such 
that z = [3w; it is clearly nonzero and positive because of positive semi-definitcness of AqLB. □ 
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We start by providing some basic Friedrichs-type energy estimates. 

Lemma A. 2. Assume u,q and Re A > solve (|A.2[) - (IA.3|1 . If e > is sufficiently small, then 
there hold the estimates 

(ReA)M^ + \ q \\ 2 + | fe || 2 < C J \U x \\u\ 2 dx (A.7) 

|ImA| J \U x \\u\ 2 dx < C J \U x \(5\u\ 2 + 8- 1 \q\ 2 )dx (A.8) 
for some C > and any 5 > 0. 

Proof. Multiply (|A.2[) by := Aq(U(x)) and take the complex L 2 product against u; taking its 
real part and denoting 

A := (A A((U(x)), L := A Q (U(x))L, 

we obtain 

(ReA)(w, A u) +Re(u, Au x ) +Re{u,Lq x ) = 0. 
Using symmetry of A and integrating by parts we get 

(Re X)(u, A u) - |Re (u, A x u) + Re (u, Lq x ) = 0. (A. 9) 

Multiply (|A.3|1 by (3 :— [3(U{x)), use (IA.6|) . take the L 2 product against q, integrate by parts and 
take its real part. This yields 

c^Wxlh +c~ 1 Mi 2 +Re(q,f3 x q) -Re{u,Lq x ) -Re(L x q,u) = 0, (A.10) 

because (3 > c _1 > 0. Since the error terms can be absorbed 

X ,L X = 0(\U X \) = 0(e 2 ), 

for e sufficiently small, and since Aq is positive definite, we obtain inequality (|A.7|) . Inequality (|A.8[) 
follows in a similar fashion, with the parameter S arising after application of Young's inequality. □ 

Corollary A. 3. There hold the estimates 

0<ReA<Ce 2 , (A.ll) 

|ImA|<Ce, (A.12) 

for some C > 0. 

Proof. Estimate (|A.11|) follows immediately from (IA.7|) . Taking 5 = e > in (1A.8|) . and using (|A.7j) 
to control |<?|^2 we can easily obtain 

(|ImA| -Ce) J \U x \\u\ 2 < 0, 

yielding (fA~T2|) . □ 
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A.l. Kawashima-type estimate. Next we carry out an energy estimate for u x of Kawashima-type 

(see nam]). 

Lemma A. 4. For each Re A > 0, A ^ 0, there holds 

KIL < C((ReA)77|u|| 2 + J \U x \\u\ 2 dx) , (A.13) 
for some C > and rj > with e 2 /ry sufficiently small. 

Proof. Denote K — K(U(x)), and take the real part of the L 2 product of Ku x against (|A.2|) . Since 
K is skew-symmetric, the result is 

Re (u x , KAu x ) = Re (X(Ku x , u» + Re (Ku x , Lq x ). (A.14) 

Noticing also that lm(Ku x ,u) = —^(K x u,u), we obtain the bound 

Re{X(Ku x ,u)) < C(ReA)(?7- 1 |M x ||2 +v\u\h) +C\lmX\ J \U x \\u\ 2 dx, (A.15) 

for any r\ > and some C > 0. We also have the estimate 

(Ku x ,Lq x ) < C{5 1 \u x \ 2 L2 +5^ 1 \q x \ 2 L2 ), (A.16) 

for any 5\ > 0, where we have used Young's inequality in both estimates. 
To estimate Re (u x , KAu x ), observe that from (|1.6[) . there holds 

Re(u x ,KAu x ) + (u x ,LBu x ) > cT 1 ^^, (A.17) 

for some c > 0. (Notice that (u x , LBu x ) € K because LB is symmetric, positive semi-definite.) 

Multiply equation (|A.3|) by L, take the L 2 product with u x and integarte by parts. This yields, 

(u x , LBu x ) = -(u xx , Lq x ) - (u x , L x q x ) - (u x , Lq). (A. 18) 

To estimate the first term, take the real part of the L 2 product of u xx against Aq times (|A.2|) . 
use A symmetric, Aq positive definite, and integrate by parts to obtain 

-Rc {u xx ,Lq x ) < —Re (X(u X7 (A ) x u)) + \{u x ,A x u x ) - Re (u x ,A x u) 

< -(Re A) Re (u x , {A ) x u) + (Im A)Im (u x , (A ) x u}+ (A.19) 

+ \ (u x ,A x u x ) -Re(u Xl A x u). 

Using (|A.7[) and (|A.8|) . and bounding the error terms (A ) X ,A X = 0(\U X \) = C(e 2 ), we get 

-Re{u xx ,Lq x ) < Ce J \U x \\u\ 2 dx + Ce\u x \ 2 L 2, (A.20) 

where the term {u x , A x u) has been bounded by 

\U x \\u\\u x \dx<^n \U x \ 3 / 2 \u\ 2 dx + J \U x \^ 2 \u x \ 2 dxj 

< Y e / ' \U x \\u\ 2 dx + je\u x \ 2 L2 . 

We also estimate 

Re(u x ,L x q x ) <Ce 2 \u x \ 2 L2 +C\q x \ 2 L2l (A.21) 

Re(u x ,Lq) <C(S 2 \u x \ 2 L2 +S^\q\ 2 L2 ), (A.22) 
for any 62 > 0, using Young's inequality. Putting all together back into (|A.18|) we get 

(u X: LBu x ) <Ce\u x \ 2 L2 +C J \U x \\u\ 2 dx, (A.23) 
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after using (|A.7[) . 

Finally, since Re A = 0(e 2 ), taking 82 — e and e 2 /r] sufficiently small, we can substitute (|A.23|) . 
(|A.15P and (|A.16[) back into (|A.17[) , absorb the small terms into the left hand side to obtain (|A.13p . 
This proves the result. □ 

Corollary A. 5. For all e > sufficiently small and Re A > 0, there holds the estimate 

{ReX)\u\ 2 L 2 + \u x \ 2 L2 <C J \U x \\u\ 2 dx 7 (A.24) 

for some C > 0. 

Proof. Take C times estimate (|A.7| and add to ()A.13|) to obtain 

C(ReX)\u\%2 + \u x \%2 < C(l + C) { \u x \\u\ 2 dx + C(Yle\)n\u\ 2 L 2. 



Take n sufficiently small, say 77 = 0(e) so that e 2 /r] remains small, and after absorbing into the left 
hand side we obtain the result. □ 

A. 2. Goodman-type estimate. Finally, we control the term J \U x \\u\ 2 by performing a weighted 
energy estimate in the spirit of Goodman [H[S] (see also 8, 25J ) . 



Lemma A. 6. Under dSOJ) - (|S2J), (|H0|) - (|H3|) . for all Re A > there holds the estimate 

Re\(\u\ 2 L3 + K|| 2 ) +C J \U x \\u\ 2 dx < Ce\u x \ 2 L2 , (A.25) 

for some C > and all e > sufficiently small. 

We first recall that there are matrices L p , R p diagonalizing matrix A such that 

(Al \ 

A := L p ARp = a p (A.26) 

\0 At) 

where A± are symmetric and positive/negative definite, and a p is scalar satisfying (|A.5[) and a p = 
0{e). Defining v := L p u, we rewrite (|A.1|) as 

Aw + Av x + Lq x = A(L p ) x RpV, 

(A.2/J 

-Qxx + q + Bv x = -B(R p ) x v, 

where 

A = L p AR p , L = L p L 7 B = BR p . 

Define 

(<t>-I P -i \ 

S := 1 (A.28) 

\ <f> + I n -p) 

where block diagonal form is in the same way as of (|A.26j) and <fr± are scalar functions of x S R 
which are bounded away from zero and satisfying 

ct>'± = Tc*\U x \ct>±, <M0) = 1 

for some sufficiently large constant c» to be determined later. Once again, we alternatively write ' 
or (-) x as derivative with respect to x. 

In what follows, we shall use (•, •) as a weighted norm defined by 

(fj) ■= (Sf,f)&. 
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With this inner product, we note that for any symmetric matrix A, 

(Af x ,f) = -~((A x + (S x /S)A)f,f) 

where S x /S should be understood as 4>'±/4>± or in each corresponding block. 
By our choice of 5* and <p±, we observe that 

({A-)' + W_/4>-)A- 
A x + (S X /S)A = a' p 

V + ^'+l<t>+)Ai / 

— c*7 



< 








(A.29) 



(Re\){v,v) + -c*(\U x \v±,v±) +-6{\U x \v p ,v p ) < -Re(Lq x ,v). (A.30) 



Proposition A. 7. Denoting v =: (v~,v p ,v+) , we obtain 

)f*(\U x \v±,v±) + - 

Proof. We take inner product in the weighted norm of the first equation of (|A.27|) against v , take 
the real part of the resulting equation, and make use of integration by parts, yielding 

(Re \)(v, v) - ((A x + (S x /S)A) x v, v) = -Re (Lq x , v) + Re {A(L' p R p v, v). (A.31) 
Noting that L' p R p ~ 0(\U X \) and the fact that A has the diagonal block (|A.26[) , we estimate 
\(AL' p R p v,v)\ < C{\U x \v±,v±) + C(\a p \\U x \v p ,v p ). 

Using this, (|A.29|) and the fact that \a p \ = 0(e) is sufficiently small and c» is sufficiently large, 
(TA~3T|) immediately yields (|A~30| . □ 

Proposition A. 8. We obtain 

(Re \){v x ,v x ) - Re (Lq x ,v) < C(0(\U x \ 2 )v, v) + n{v x ,v x ) (A.32) 
for sufficiently small n > 0. 

Proof. We now take the inner product of the derivative of the first equation of (|A.27[) against v x ■ 
We thus obtain 

M v x,v x ) + ((Av x ) x ,v x ) + {(Lq x ) x ,v x ) = ((AL' p R p v) x ,v x ) (A. 33) 

where we estimate by integration by parts, 

((Av x ) x ,v x ) = (A x v x ,v x ) - ^((A x + (S x /S)A)v x ,v x ) = (0(\U x \)v x ,v x ) 
((AL' p R p v) x ,v x ) = (0(\U x \)v x ,v x ) + (0(\U x \)v,v x ) 
and by using the second equation and the semi-definite condition LB > 0, 
((Lq x ) x ,v x ) = (Lq xx ,v x ) + (L 

= (L(q + Bv x + BR' p v),v x ) + (L x q x ,v x ) 

= -(Lq x ,v) - ((L x + (S x /S)L)q,v) + (LBv x ,v x ) + (LBR' p v,v x ) + (L x q x ,v x ) 
> -(Lq x ,v) + (LBR' p v,v x ) + (L x q x ,v x ) - ((L x + (S x / S)L)q,v). 
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Thus, fQ3l) yields 

(Re A) (v x , v x ) - Re (Lq x , v) 

_ (A.34) 

< (0(\U X \ )v,v) + n(v x ,v x ) + (L x q x ,v x ) ~ ((L x + (S x /S)L)q,v). 

By testing the second equation against q, it is easy to see that 

) + (q,q) < C(v x ,v x ). 

Thus, we have 

(L x q x ,v x ) - ((L x + (S x /S)L)q,v) < C(v x , v x ) 1 ' 2 {{0{\U x \ 2 )v x , v x ) + (0(\U x \ 2 )v, «)) 1 * 

Using the standard Young's inequality and absorbing all necessary terms into the right hand side of 
(|A.34[) . we thus obtain from (|A.34|1 the important estimate, (|A.32[) , which proves the proposition. □ 

Combining Propositions I A. 7l and [A~8l we are now ready to give: 

Proof of Lemma [A . 61 Adding (|A.32[) with (|A.30[) , noting that the "bad" term Re (Lq x , v) gets can- 
celed out, and using the fact that \U X \ = 0(e) is sufficiently small, we easily obtain 

ReX({v,v) + (v x ,v x ))+6(\U x \v,v) < r/{v x ,v x ) (A.35) 
which by changing v to the original coordinate u yields the lemma. □ 
Proof of Theorem [TT6l Add Ce times (|A~24| to ((A~25l) to get 



(ReA)(l + Ce)\u\i 2 + {C + CCe) J \U x \\u\ 2 dx < 0, 
which readily implies Re A < 0, yielding the result. □ 



Remark A. 9. Theorem If .61 can be extended to the non-convex case, that is, when the principal 
characteristic mode is no longer genuinely nonlinear (hypothesis (|H2j) does not hold). For that 
purpose, it is possible to modify the Goodman-type weighted energy estimate by means of the 
Matsumura-Nishihara weight function w [28] (introduced to compensate for the loss of monotonicity) , 
satisfying 

-\(wa p + w x ) = \U X \, 

which replaces the 1 in the weight matrix function S in (|A.28|) . This procedure was carried out for 
the viscous systems case by Fries [2] and it can be done in the present case as well at the expense 
of further book-keeping. Note that the existence result of [TTJ, [TS] includes non-convex systems, a 
feature that might be useful in applications. 

Appendix B. Pointwise reduction lemma 

Let us consider the situation of a system of equations of form 

W x =A e {x,\)W, (B.l) 

for which the coefficient A £ does not exhibit uniform exponential decay to its asymptotic limits, but 
instead is slowly varying (uniformly on a e-neighborhood V, being e > a parameter). This case 
occurs in different contexts for rescaled equations, such as (|2.f 6[) in the present analysis. 

In this situation, it frequently occurs that not only A e but also certain of its invariant eigenspaces 
are slowly varying with x, i.e., there exist matrices 



L e - (*)> K e = (Rl #2) (*) 
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for which L e M e (x) = / and |LM'| = |L'R| < CS e (x), uniformly in e, where the pointwise error bound 
<5 e = S e (x) is small, relative to 

M* :=UA<R'(x) = ( M J ^(x) (B.2) 

and ""' as usual denotes d/dx. In this case, making the change of coordinates W e = WZ, we may 
reduce (|B.1|) to the approximately block-diagonal equation 

Z e/ = M c Z e + <5 e e e Z e , (B.3) 

where M e is as in (|B.2[) . Q e (x) is a uniformly bounded matrix, and 8 e {x) is (relatively) small. 
Assume that such a procedure has been successfully carried out, and, moreover, that there exists an 
approximate uniform spectral gap in numerical range, in the strong sense that 

min <j(Rc A/j 5 ) — max cr (Re M|) > rf( x )i f° r & U x > 

with pointwise gap r/^^x) > t]q > uniformly bounded in x and in e; here and elsewhere KeN := 
7}(N + N*) denotes the "real", or symmetric part of an operator N. Then, there holds the following 
pointwise reduction lemma, a refinement of the reduction lemma of [23] (see the related "tracking 
lemma" given in varying degrees of generality in [3l 122 ] I3T1 [41] [38] ) . 



Proposition B.l. Consider a system (|B.3[) under the gap assumption |B|) ; with e uniformly 
bounded in e S V and for all x. If, for all 6 6 V, sup a . 6R (<5 e /ry e ) is sufficiently small (i.e., the ratio of 
pointwise gap rj e {x) and pointwise error bound 6 f ~(x) is uniformly small), then there exist (unique) 
linear transformations $| (x, A) and ^(^j ty> possessing the same regularity with respect to the var- 
ious parameters e, x, A as do coefficients M e and S e (x)Q t (x), for which the graphs {(Zi,$|(Zi))} 
and {(<fr\(Z2), Z2)} are invariant under the flow of (|B.3|) . and satisfying 

sup I & I < Csup(57ry c ). 

R R 

Moreover, we have the pointwise bounds 

m(x)\ <C [ X e-^^ {z)dz 5^y)dy, (B.4) 

and symmetrically for $f . 



Proof. By a change of independent coordinates, we may arrange that r) e (x) = constant, whereupon 
the first assertion reduces to the conclusion of the tracking/reduction lemma of [23] . Recall that 
this conclusion was obtained by seeking $| as the solution of a fixed-point equation 




^6%y)Q(^ 2 )(y)dy. 



Observe that in the present context we have allowed d e to vary with x, but otherwise follow the proof 
of [53] word for word to obtain the conclusion (see Appendix C of [23 , proof of Proposition 3.9). 
Here, Q($|) = 0(1 + |$|| 2 ) by construction, and \T V ^ X \ < Ce^^-vi Thus, using only the fact 
that I $|| is bounded, we obtain the bound (]B.4|) as claimed, in the new coordinates for which rf is 
constant. Switching back to the old coordinates, we have instead \T^ X \ < Ce^ ^ v ' (z)dz , yielding 
the result in the general case. □ 
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Remark B.2. From Proposition lB.il we obtain reduced flows 

Zf = M{Z{ + 6 e (Q n + ®y&l)Zl 

z% = m\z\ + «5 £ (e 22 + eii$f)z|. 

on the two invariant manifolds described. 
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